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Abstract 

The signal generated by a weak impulse propagates in an oscillatory way and 

dispersively in a gravitationally compacted granular chain. For the power-law 

type contact force, we show analytically that the type of dispersion follows 

power-laws in depth. The power-law for grain displacement signal is given 

by h 4V p' where h and p denote depth and the exponent of contact force, 

and the power-law for the grain velocity is h 4X3 p. Other depth-dependent 

power-laws for oscillation frequency, wavelength, and period are given by com- 

bining above two and the phase velocity power-law /i 2 p. We verify above 

power-laws by comparing with the data obtained by numerical simulations. 
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Physics of granular materials attracts great interest recently [|]] , since these materials are 
ubiquitous around us and their properties are unique and also useful in many applications 
|2|j3|] . The propagation of a sound or a weak elastic wave in a granular medium is also one of 
interesting subjects related to the properties of granular matter A rather simple system, 
the granular chain with Hertzian contact ]5J, has been revived by finding a soliton in trans- 
mitting elastic impulse. This soliton, existing in a highly nonlinear regime of a horizontal 
Hertzian chain was first predicted by Nesterenko |J and its experimental verification was 
performed by Lazaridi and Nesterenko [|7| and recently by Coste et al. ||. Even though 
three-dimensional granular systems may not follow simple Hertzian contact force law due to 
geometrical effect ||, the one-dimensional granular chain with nonlinear contact force is still 
interesting. It may describe a fundamental feature of the dynamics of nonlinear granular 
chain which appears in many areas of nature. In addition, one- dimensional system is usually 
the starting point of studying higher dimensional systems. 

It is well-known that the velocity of an elastic impulse scales as P 1 / 6 or h 1 ^ 6 for the 
Hertzian chain [EJ , where P is the pressure, linearly proportional to the depth h for vertical 



chain. Sinkovits and Sen [TO] extended this to arbitrary nonlinear contact force of power-law 

type F oc 5 P , where 5 denotes overlapped distance between adjacent grains. They showed 

fi— -)- 

that the signal velocity v p h scales as /r p' 2 for p > 1 at large h. This has been simply 
obtained by considering the well-known relation v p h oc y/Jl, where fi is the elastic constant 
which is given by fi ~ h l ~ p for the above power-law type contact force. As far as we know, 
however, no power-law dependences on depth of the signal characteristics, such as oscillation 
frequency, period, and wavelength have been found in the gravitationally compacted chain. 

In this work, we study the propagation of acoustic or weak impulses in the gravitationally 
compacted granular chain. We derive analytically the power-law behaviors of signal charac- 
teristics which depend on depth or time. We treat here a rather weak impulse which makes 
grain motion oscillatory and can be treated analytically even though it contains nonlinearity. 
The other extreme which is a highly nonlinear regime has been studied by Nesterenko ||. 
Initial impulse may be used as a parameter which controls the solitariness of signal. The 



power-law behaviors for a wide range of impulse will be discussed in a separate work JTT . 

We would like to obtain analytically the exponents of various power-laws, such as grain 
displacement, grain velocity, and oscillation period, frequency, and wavelength. We first 
solve the equation of motion of a grain displacement under gravity in the small oscillation 
or weak impulse regime in which the equation of motion under gravity can be mapped into 
the equation for the horizontal linear chain with varying force constant at each contact. 
The normal mode solution of the equation of motion can be obtained analytically in the 
continuum or long wavelength limit. The asymptotic behavior of the normal mode gives rise 
to the correct power-law behavior in depth, since the equation of motion has been changed 
into a linear form. Once we get the information on the grain velocity, all sorts of power-laws 
mentioned above can be obtained. Since the equation of motion for grain velocity is not 
linear, the normal mode solution may not work to obtain power-law in depth. Therefore, 
we construct fully nonlinear forms describing displacement and velocity signal and obtain 
their depth-dependence behaviors. Our solution is quite general and gives rise to generic 
power-laws for arbitrary exponent p of the contact force in the oscillating regime. 

The equation of motion of n-th grain at z n is given by 

mz n = 7][{A - (z n - z n _i)} p - {A - (z n+ i - z n )Y] 

+ mg, (1) 

where z n is the distance from the top of chain to the center of the n-th spherical grain, m is 
the mass of grain, Ao is the distance between adjacent centers of the spherical grain, and r\ is 
the elastic constant of grain. Therefore, the overlap between adjacent grains at nth contact 
is 5 n = A — (z n+ i — z n ). It is usually impossible to solve general nonlinear problems in an 
analytical way. Therefore we may not solve the nonlinear differential equation of Eq. (1) 
exactly. But we may treat it analytically in a small oscillation regime which can be achieved 
by applying a weak impulse. 

For this purpose, we introduce a new variable 
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n (mal\ 1,p 

ip n = z n - nA + J2 [~yJ > (2) 

where the last term is the sum of overlaps up to the nth grain and set z = ip = 0. This 
change of variable makes Eq. (1) into an equation for the linearized horizontal chain with 
varying force constant at each contact, i.e., 

d 2 

m—1p n = -fl n (lp n - V'n-l) + U n+ iOn+l - Ipn) (3) 

where « n = uin 1- ? is the force constant of the nth contact and ui = mpg V is the 

force constant of the first contact. Use of the condition of small oscillation 

|V>n - ^n-l| < f^pl ( 4 ) 

has been made to obtain Eq. (|3]) approximately. The expression of Eq. (|3|) in the continuum 
limit, i.e. the lattice constant a = 5h — > 0, is given by 



P 9 l ft, A i -l 9 



(5) 



where fj,(h) = h l p denotes the depth-dependence of force constant, and p = m/a and 
7~i = \i\<a are the linear density and the tension of a chain at the first contact, respectively. 



We set c\ = Jtx/p which is the well-known speed of wave in the string of tension n and 
line density p. 

i _ i 

We now choose ip^(h,t) = u^(h)e~ tCia v ^ as a normal mode solution. Then the depth- 
dependent function u^(h) satisfies 



which is a type of Bessel's differential equations [12]. If we consider a solution propagating 



to the positive /i-direction, the solution of Eq. (|6|) is given by the Hankel function |T2 



u ( (h) = htHU(6W), (7) 



w here£ = ^,7 = f + £=§(l + £),0 



The asymptotic form of Eq. fl7|) at large h for a fixed v is 



Uc (h) w ^/A^-ie^^-l-f] ( 8 ) 
and the displacement function becomes 

r 1 - 1 i 

^ c (/i,t) ~/i f_ ie L 1 J (9) 

The amplitude of displacement signal is given by the envelope function of the asymptotic 
solution, which scales as 

A(h) ~ h*-l = /T^ 1 "*) (10) 
The phase velocity v p h is obtained by setting the phase of Eq. @ constant, i.e., ^h 1 — 

_L 1 

Cia 2 p 2 (t =constant. We obtain 

^ = | =Cl0 i«-), i M'-y (id 

The group velocity scales as that of phase velocity and the dispersion relation is linear but 
depth-dependent in this case. 

Since Eq. (||) is a linear differential equation, the envelope function of a normal mode 
solution in Eq. ([|) can describe the depth-dependence of the displacement signal which may 
be given by the linear combination of all normal modes of different frequency. Therefore, 
one can expect that the depth-dependent behavior predicted by Eq. ( |10"1) may agree well 
with the data given by numerical simulation which will be shown in what follows. 

The normal mode solution of Eq. @ alone, however, cannot give appropriate predictions 
on the changes in frequency, wavelength, and period of the signal as it propagates down. 
The information on the grain velocity may give the characteristics of signal dispersion by 
combining it with the displacement and velocity of signal. To obtain the depth-dependence 
of grain velocity we write Eq. (I5D as 



9 n ^ d 



d 



(12) 
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We set the constant factor of Eq. (|5]) unity for our convenience, since it has nothing to do 
with depth-dependence behavior. To draw the depth-dependent behaviors of grain velocity, 
frequency, etc. we set the frequency uj(h) ~ h a , the displacement function 

V>(M) ~ hr^-v) e ik{h)h ~ ilAj{ - h) \ (13) 

and the grain velocity 

v{h,t)^h p e lk{h)h - lw{h)t+ ^, (14) 

where a and (3 will be determined and cb is the phase difference between displacement and 
velocity signal. The depth-dependent wavenumber k(h) is given by k(h) = v p h{h)/uj(h). 

Since we have two unknowns a and (3 to be determined, we need two independent equa- 
tions for these. One is given by Eq. (|T^) and the other the relation u(h) ~ v(h)/ip(h). 
Substituting Eqs. (fL3|) and (|b|) into Eq. ( ^2|) gives rise to (3 = ~ — ^ — 2a and the relation 
uj{h) ~ v{h)/ip{h) yields a = P + \ — We obtain the power-law behaviors of grain velocity 
and frequency from these two equations as follows: 

v(h) oc /T^ + p) (15) 
w(h) oc hJ~^ (16) 

The characteristic time of oscillation which is expressed by the period is given by the 
inverse of frequency or the ratio of displacement to grain velocity, i.e., 

T(h) = ^ oc u(h)- 1 ~ (17) 

The characteristic length of oscillation, on the other hand, which is expressed by wavelength 
is given by multiplying T(h) by phase velocity, i.e., 

\{h) = T(h)v ph (h) ~ h 1 / 3 (18) 

Interestingly enough, characteristic length of oscillation does not depend on contact force 
within the linear approximation of Eq. (3). 



We now compare above results with molecular dynamics simulations performed for Eq. 
(1) for arbitrary p. To perform numerical simulation for Eq. (1), we choose a vertical chain 
of N = 2 x 10 3 grains and neglect plastic deformation. As a calculational tool, we use the 
third-order Gear predictor-corrector algorithm ||13|| . We choose 10 _5 m, 2.36 x 10 _5 kg, and 
1.0102 x 10~ 3 s as the units of distance, mass, and time, respectively. These units gives the 
gravitational acceleration g — 1. We set the grain diameter 100, mass 1, and the elastic 
constant r\ of Eq. (1) is given by rj = (1 + p)b, where b depending on modulus is chosen as 
5657 for this molecular dynamics simulation. The equilibrium condition 

mgn = rjS^ (19) 

has been used for the (n + l)th grain of a vertical chain. Even though a criterion |8|,|l4j| for 
initial impulse to neglect plastic deformation and viscoelastic dissipation in experimental 
situation, we do not care about that criterion for the numerical simulation. For the purpose 
of this work, however, we choose a rather weak initial impulse v j = 0.1 in our program units. 
There is a regime of initial impulse in which the signal follows the same power-laws. This 
will be shown in a separate paper JTT] . 



Figure 1 shows the snap shots of amplitudes (a) and corresponding grain velocity signals 
(b) propagating down the vertical chain with Hertzian contact (p = 3/2). The leading 
amplitude of displacement of each signal in Fig. 1(a) corresponds to the leading part of each 
velocity signal in Fig. 1(b). 

We focus on the leading amplitudes of displacement and velocity signal and plot them 
in logio-logio scale in Fig. 2 which shows that both displacement and velocity peak de- 
creases in power-laws of depth. The explicit expressions for the depth-dependent behaviors 
of leading amplitudes of displacement and velocity are given by A max (h) oc Z;,- - 083 ^ - 0003 
and v max (h) oc /;,- - 2500±00001 . We also obtain other depth-dependent power-laws showing 
dispersiveness of the signal. They are the elasped time to reach to maximum amplitude, 
T ma x(h), which describes the period and the number of particles participating at the leading 
part of velocity signal, N(h), which describes the wavelength. The power-law exponents of 



these quantities with error bounds are T max ~ / l o.i7o±o.oo2^ ^ ^o.338±o.oo4_ ^^ggg values 
are in good agreement with theoretical predictions for the Hertzian p = 3/2. 

We obtain peak values of displacement and grain velocity signal for other values of p 
and plot them in Fig. 3. One can see a very nice fit to the theoretical curves up to p = 2. 
For large values of p, the deviation from theory occurs especially in grain velocity. This 
is understandable because nonlinearity becomes stronger as p increases and grain velocity 
contains more nonlinearity than displacement. 

In conclusion, the propagating feature in vertical chain is dispersive due to gravity even 
though total energy and momentum are conserved. The effect of gravity induces the change 
in force constant at every contact. Therefore, the signal is no longer a soliton which is the 
propagating mode in horizontal chain |||§. We treat the problem analytically for arbitrary 
power-law type of nonlinear contact forces and obtain general features of dispersive phenom- 
ena for a weak impulse in a gravitationally compacted chain. The normal mode solution for 
displacement has been obtained in the small oscillation and continuum limit. This normal 
mode solution describes the depth-dependent power-law of displacement signal. We find 
various power-laws describing signal characteristics depending on depth. This dispersive 
property is obtained by constructing both displacement and grain velocity function appro- 
priately. We also perform numerical simulations for various power-law type contact force 
for a rather weak impulse v j = 0.1 and show that the results agree well with our theoretical 
predictions. 

The general features of soliton damping due to gravity may be given by studying similar 
work for a wide range of impulse, which will be given in a separate paper jTT|. The properties 



of signal propagation studied in this work are fundamentals of the dynamics of granular 
chain under gravity. This work may be extended to higher dimensions and to more practical 
models for applications. 
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Figure Captions 

Fig. 1 : (a) Snap shots of displacement of the propagating wave in a gravitationally com- 
pacted granular chain with Hertzian contact force law. Initial impulse is Vi — 0.1. (b) 
Snap shots of grain velocity corresponding to (a). 

Fig. 2 : Logio-logio plots of leading peaks of displacement and grain velocity shown in 
Fig. 1. Slopes of the straight lines are —0.0835 (displacement) and —0.250 (velocity), 
respectively. 

Fig. 3 : Comparison of theory and simulation for power-law exponents for various p. Sim- 
ulation data are obtained for the leading peaks of displacement and grain velocity 
signal. 
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